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Abstract. This is a continuation of our earlier study of the shape parameter c contained in the 
famous multiquadrics (— 1)^1 (c 2 + ||a;|| 2 ) /3 , /3 > 0, and the inverse multiquadrics (c 2 + ||a;|| 2 ) /3 , /3 < 0. 
In [5] and [TU] we present criteria for the optimal choice of c, based on the exponential-type error 
bound [BJ. In this paper a new set of criteria is developed, based on the improved exponential- type 
error bound[?J. Unlike [UJ and [TU], where the range of c is only a subset of (0, oo), the range of c 
we investigate here is the entire interval (0, oo). The drawback is that the distribution of the data 
points is not purely scattered. However it seems to be harmless. 
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1 Introduction 

We begin with some basic ingredients of our theoretical ground. 

Let T n denote the n-simplex in R n whose definition can be found in [3J. A 1-simplex is a line 
segment, a 2-simplex is a triangle, and a 3-simplex is a tetrahedron with four vertices. 

Let Vi, l<i<n+lbe the vertices of T n . Then any point x € T n can be written as convex 
combination of the vertices: 

n+l n+l 

X = ^2 °i v i> X! Ci = 1 ' Ci - °' 

z=l i=l 

The numbers ci, • • • , c n +i are called the barycentric coordinates of x. 

For any n-simplex T n , the evenly spaced points of degree I are those points whose barycentric 
coordinates are of the form 

(-J-, ~f~'" ' ' )' nonnegative integers with k% = I. 

If we let P™ denote the space of polynomials of degree not exceeding I in n variables, it's easily seen 

that the number of evenly spaced points of degree I is exactly N = dimP" = f n ^ ' V Also, such 

points form a determining set for P™, by [2]. 

In this paper the interpolation will happen in n-simplex and the set X of centers(interpolation 
points) will be the evenly spaced points in the n-simplex. 
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The radial function we use is 

h(x) := r(-^)(c 2 + \\x\\ 2 )?, (i G i?\2iV> , c> (1) 

, where ||x|| is the Euclidean norm of a; € R n , T is the classical gamma function, and j3, c are 
constants. Note that this definition is slightly different from the one mentioned in the abstract. We 
adopt (1) because it will greatly simplify its Fourier transform and our future work. The function 
h{x) in (1) is conditionally positive defmite(c.p.d.) of order m = max{0, [§]} where ["£] means the 
smallest integer greater than or equal to 4. Further details can be found in jllj . 

Given data points (xj,yj), j = 1, • • • ,N, where X — {x\, • • • , xn} is a subset of R n and z/j's are 
real or complex numbers, our interpolant will be of the form 

N 

s(x) = p(x) +2_, c jh{x — Xj) (2) 
3=1 

, where p{x) is a polynomial in P^_ 1 to be determined and Cj's are coefficients to be chosen. 

As is well known in the theory of radial basis functions, if X is a determining set for P£ l _ 1 , there 
exists a unique polynomial p(x) and unique constants ci, • • • , cjv satisfying the linear system 

N 

p(xj) + >^ Cjhjxj - Xj) = yi, i = l,---,N 

3=1 

(3) 

N 

^c 3 q{x 3 ) = 

3=1 

where q ranges over all basis elements of P^-i- All these can be found in jTTj . 



1.1 Fundamental theory- 
Each function h of the form (f) induces a function space Ch, m , called native space, whose definition 
and characterization can be found in [TT] . [1 2j . [4] . [5] . [8] . and [T4] . Here m — max{0, [4 ]}. Also, there 
is a seminorm \\f\\h for each / <E Ch. m - In our theory every interpolated function belongs to the 
native space. 

Before entering the main theorem, let's introduce two constants. 
Definition 1.1 Let n and f3 be as in (1). The numbers p and Aq are defined as follows. 

(a) Suppose < n — 3. Let s = |" "~f~ 3 ] . Then 

(i) if? < 0, p= and A = (glflMil^ ; 
(ii) ifP>0, p=l + ^ and A = ^+*+>KW+>h<*>»+*) 
where m = [^] . 

(b ) Suppose n — 3</3<n — 1. Then p = 1 and Aq = 1 . 



(c) Suppose P>n-1. Let s = - \^^\ . Th 



en 



p=l and An = -; n ; ; r where m = \— 1 . 

H (2m + 2)(2m+l)---(2m-s + 3) 1 2 1 
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Our criteria for the optimal choice of c is based on the following theorem which we cite directly from 
[7] but with a slight modification to make it easier to understand. 

Theorem 1.2 Let h be as in (1). For any positive number bo, let C = max{ g 2 ^-, 8p} and 5q = 
For any n-simplex Q of diameter r satisfying ^ < r < (note that ^ < bo), if f € Ch. m , 

\f(x) - s(x)\ < 2^7T^ 1 y /Ec7; i c^ l ^V3CVS(X / )i\\f\\ h (4) 

holds for all x € Q and < 5 < Sq, where s(x) is defined as in (2) with x±, • • • , xn the evenly spaced 
points of degree I in Q satisfying -^jg < I < . The constant a n denotes the volume of the unit 
ball in R n , and < A' < 1 is given by 




which only in some cases mildly depends on the dimension n. 

Remark: This seemingly complicated theorem is in fact not difficult to understand. Note that the 
right-hand side of (4) approaches zero as S tends to zero. Also, (4) is greatly influenced by the shape 
parameter c. However, in order to make it useful in the choice of c, we still have to do some work. 
We begin with the following definition. 

Definition 1.3 For any a > 0, the class of band-limited functions f in L 2 (R n ) is defined by 

B a = {fe L 2 (R n ) : /(£) - if |£| > a} 

, where f denotes the Fourier transform of f. 
Now we cite Theoreml.6 of [9] as a lemma. 

Lemma 1.4 Let h be as in (1) with j3 > 0. Any function f in B a belongs to Ch, m and 

H/IU < y/m\S(m,n)2- n -^K- n -ia^e^c^\\f\\ La{Rn) 
where c,f3 are as in (1) and S(m,n) is a constant determined by m and n. 

Corollary 1.5 Let h be as in (1) with ft > 0. If a > and f £ B a , the inequality (4) can be 
transformed into 

\f(x) - s(x)\ < 2 - 2 -f" 7r "^V^V^^V / ™!5(m,n)a^c i± ^-'e^^(A')^|/||L 2 (fl") (5) 

In orde to handle the case j3 < 0, we need the following lemma which is just Theoreml.7 of [5]. 

Lemma 1.6 Let h be as in (1) with j3 < such that n + /3 > 1 or n + /? = — 1. Any function f in 
B a belongs to Ch m and satisfies 

ll/IU < 2- n - — TT- n —*a^e-c^\\f\\ L 2 {Rn) . 

Corollary 1.7 Let h be as in (1) with f3 < such that n + f3 > 1 or n + f3 = —1. If a > and 
f G B a , the inequality (4) can be transformed into 

\f(x) - s(x)\ < 2- 2 -^V^V^VA^V3C^ i±4±i c i± ^- z e^^(A')'||/|U 2 (^) (6) 
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Note that Corollary 1.7 does not cover the frequently seen case (3 = — l,n= 1. For this we need the 
following lemma which is just Lemma2.1 of [S]. 

Lemma 1.8 Let h be as in (1) with f3 = — l,n = 1. For any a > 0, if f € B a , then f G C/i, TO cind 

( } 1/2 

ll/H, < (27r)-"2-i -i- / |/(0| 2 d£ + - / |/(e)| 2 V^Te c l«l^ 
^ /C (l) J| € |<i a Ji<| 5 |< CT J 

if - < a, where an = —K=, and 

J c ' u 2V3 

»/»^«-" 2 -'{icSl)l,J /({)|2rfS } ,/2 

i/ i > a-. 

Corollary 1.9 Lei h be as in (1) with (3 — — l,n = 1. If a > and f € B a , the inequality (4) can 
be transformed into 

\f{x) - s(x)\ < 2- 2+ ^^n^^^n^^V3Cc?- l VS(\ , )*{A + B} 1 ' 2 (7) 

where A = ^ \f(£)\ 2 d£ for all c, B = 2^1^^ |/(0| 2 V^ eC|?l ^ if \ < v, and B = Q 

if l >o-. 

J c — 



2 Criteria of Choosing c 

The results of Sectionl provide us with useful theoretical ground for choosing c. Note that in the 
right-hand side of (5), (6) and (7), there is always a main function determined by c. Let's call it MN 
function, denoted by MN(c), as in [ID] . Its graph is called MN curve. Then finding the optimal 
value of c is equivalent to finding the minimum of MN(c). Let's analyze as follows. Here the range 
of c is the entire interval (0,oo). 



Casel. f3 > and n > 1 Let / € B a and h be defined as in (1) with /3 > and n > 1. Un- 
der the conditions of Theorem!. 2, for any fixed S satisfying < S < 5q, the optimal choice of c is as 
follows. 

(a) If 1 + (3 — n — 41 > 0, where / is defined as in Theorem!. 2, let c be as small as possible. 

(b) If I + (3 — n — 41 < 0, where / is defined as in Theoreml.2, the optimal c is where the 
minimum of MN{c) :— e~c~ ' occurs. 

Reason: (a)By (5) the MN function is MN(c) = e^c ^ 3 ^"' 4 ' . If 1 + ft - n - 41 > 0, MN(c) is an 
increasing function. Hence the smaller c is, the better it is. (b)If 1 + f3 — n — 41 < 0, MN(c) —> oo 
both as c — ¥ + and c — >■ oo. The minimum of MN(c) corresponds to the optimal choice of c because 
it makes the error bound (5) minimum. 

Remark: For (b), the optimal c can be easily found by Mathematica or Matlab. 



Numerical Results: 
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Graph of the MN Curve with (5=0.01 and 1=5 

MN(c) 



0.025 




0.020 




0.015 




0.810 






10 20 


Figure 1 : Here n = 1, /3 = 1, bo = 1 and a = 1 . 
Graph of the MN Curve with <5=0.008 and 1=6 



MN(c) 
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Figure 2: Here n = 1, /3 = 1, &o = 1 and cr = 1. 
Graph of the MN Curve with (5=0.006 and 1=7 
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Figure 3: Here n = 1, /3 = 1, 6q = 1 and a = 1. 
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Graph of the MN Curve with (5=0.004 and 1= 1 1 
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Figure 4: Here n = 1, j3 = 1, bo = 1 and er = 1. 
Graph of the MN Curve with (5=0.002 and 1=21 

MN(c) 
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Figure 5: Here n = 1, /3 = 1, &o = 1 and a = 1. 
For /3 < 0, we separate it into two cases. 



Case2. 



/3 < and n + j3>lovn + j3 



-1 



Let f £ B a and /i be denned as in (1) with /3 < and 
n + f3 > 1, or n + /3 = —1. Under the conditions of Theoreml.2, for any fixed 5 satisfying < 5 < Sq, 
the optimal 
Theoreml.2. 



the optimal value of c is the number minimizing MN(c) :— e^c + ^ ■* where / is defined as in 



Reason: This is an immediate result of Corollaryl.7. Note that in this case 1 + (3 — n — 4/ is 
always negative. So MN(c) — > oo both as c — > + and c — > oo. 

Remark: Case2 covers (3 = — 1 and n > 2. 

Numerical Results: 
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Graph of the MN Curve with 6 
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Figure 6: Here n = 2, (3 = — 1, &o = 1 and cr = 1. 
Graph of the MN Curve with 5=0.008 and 1=6 

MN(c) 




Figure 8: Here n = 2, /3 = — 1, &o = 1 and a = 1. 
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Graph of the MN Curve with <5=0.004 and 1=1 1 

MN(c) 
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Figure 9: Here n = 2, j3 = —1, bo = 1 and a = 1. 
Graph of the MN Curve with 5=0.002 and 1=21 

MN(c) 




Figure 10: Here n — 2, /3 = — 1, £>o = 1 and a = 1. 



Now we begin the case /3 



-1 and n = 1. 



Case3. 







-1 and n 



Let / G Bo- and /i be defined as in (1) with ft = — 1 and n = 1. 
Under the conditions of Theoreml.2, for any fixed 5 satisfying < 5 < So, the optimal value of c is 
the number minimizing MN(c), defined by 



MiV(c) 



C2 



K (l) 



2^3- 



if 0<c< i, 
if i < c 



, where I is defined as in Theoreml.2 and ICo is the modified Bcsscl function. 

Reason: In (7), the part of the upper bound influenced by c is c%~ l {A + B}i . Note that for 
< c < ~, we have a < - and 



.4 



£o(l) 7|{|<i 



\m)\ 2 dc 



Mi) 



l/(0l 2 <*e = F7ri 



l/lli 2 (fl") 
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because / € B a . Therefore 

J- l {A + B}-> = c^-'VI = - J_ c #-'||/|| £ , (Jp .j 

if < c < ±. Now, if i < c, 

^o(i) yifi<i ii<ifi<<T 



< C 2 ' 



^o(i)7| 4 |<i A<l£l<" 

~| 1/2 

^oUjJ|4i«7 ^iei<<7 I 



r i i 1 / 2 



= c 2 



{oi) +2 ^ e 1 11/11 w 



As a result, the MN function is 

1 ifO<c<i, 



MJV(c) 



c4 "'{d(T)+ 2 ^^ eCC7 }' if i< c - 



Numerical Results: 
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Graph of the MN Curve with (5=0.008 and 1=6 

MN(c) 
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Figure 12: Here n = 1, j3 = — 1, &o = 1 and £7 = 1. 
Graph of the MN Curve with c5=0.006 and 1=7 
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Figure 13: Here n = 1, /3 = — 1, &o = 1 and £7 = 1. 
Graph of the MN Curve with <5=0.004 and 1=1 1 
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Figure 14: Here n = 1, f$ = — 1, fro = 1 and a = 1. 
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Graph of the MN Curve with <5=0.002 and 1=21 
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Figure 15: Here n = 1, /3 = —1, &o = 1 and a = 1. 

Note that the optimal c increases rapidly as i5 becomes small. 

Final Remark: Both [9] and this paper deal with the interpolation of band-limited functions. 
In [9] the range of c is [12pY / ne 2 ™ 7 "7„(m + 1)6, oo) and is only a subset of (0,oo). In order to 
make the left endpoint of the closed-open interval small enough, often S must be very small. The 
consequence is that a huge number of data points will be involved, making the criteria of choosing c 
only theoretically valuable, especially for n > 2. Now we have thoroughly removed this restriction. 
In this paper the range of c is (0, oo). 
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